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In whispering gallery mode (WGM) resonators light is guided by continuous total internal re- 
flection along a curved surface. Fabricating such resonators from an optically nonlinear material 
one takes advantage of their exceptionally high quality factors and small mode volumes to achieve 
extremely efficient optical frequency conversion. Our analysis of the phase matching conditions 
for optical parametric down conversion (PDC) in a spherical WGM resonator shows their direct 
relation to the sum rules for photons' angular momenta and predicts a very low parametric os- 
cillations threshold. We realized such an optical parametric oscillator (OPO) based on naturally 
phase-matched PDC in Lithium Niobate. We demonstrated a single-mode, strongly non-degenerate 
OPO with a threshold of 6.7 fiW and linewidth under 10 MHz. This work demonstrates the remark- 
able capabilities of WGM-based OPOs and opens the perspectives for their applications in quantum 
and nonlinear optics, particularly for the generation of squeezed light. 

PACS numbers: 42.65.Yj, 42.60.Da 



WGM resonators have found applications in many ar- 
eas of photonics, including quantum and nonlinear op- 
tics, spectroscopy, biophysics and optomechanics (see [l|- 
Q and references therein). These areas benefit from a 
high quality factor (Q) and strong field confinement of 
WGM resonators in two different aspects. First, WGM 
resonators provide strong coupling of the confined field 
to other systems such as quantum dots 0, atoms 
mechanical oscillators [6] or other cavities Q. Second, 
WGM resonators provide strong nonlinear coupling be- 
tween optical fields. Here the possibilities lie in utilizing 
either third order nonlinear effects in fused silica [8|4lQ|, 
or nonlinear effects in crystalline materials [lll-Hij. 

In this letter we demonstrate optical parametric os- 
cillations with extremely low thresholds and narrow 
linewidths in a crystalline WGM disk resonator. This is 
achieved by taking advantage of the strong second-order 
optical nonlinearity of Lithium Niobate and of the high 
quality factor of our WGM resonator, limited by the ma- 
terial absorption only. This research is motivated by the 
perspective to develop a miniature, stable, narrow-line 
tunable OPO for spectroscopy applications, to study the 
OPO dynamics far above the threshold, and by the wide 
and successful use of conventional OPOs for the genera- 
tion of non-classical light. 

To achieve efficient parametric processes in a nonlin- 
ear material phase matching conditions have to be ful- 
filled. For parametrically interacting plane waves in a 
bulk material, these phase matching conditions corre- 
spond to conservation of energy and linear momentum 
of pump (p), signal (s) and idler (i) photons: 

huj p = huo s + fthJi, hk p = hk s + Tiki. (1) 
In spherical geometry the orbital angular momentum of 



the photons is conserved, along with the energy. This 
gives rise to WGM selection rules for the PDC pro- 
cess and determines the OPO threshold in-coupled pump 
power for various modes [lEj : 



Pth = 



XpQp \167rcrx^Qs 



(2) 



where c is the speed of light, X p is the pump wavelength 
n = n e (X p ) = n (X s ) is the phase-matched refractive 
index, and is the effective quadratic nonlinearity. 
Let us point out that the OPO threshold power (j2j) is 
closely related to the second harmonic saturation power 
Po = 4Pth P3, provided that the same WGMs are used. 

The factor a in (0) is determined by the overlap of 
the pump, signal and idler WGM normalized eigenfunc- 
tions ^p, s ,i- I n spherical geometry these functions are ex- 
pressed via spherical Bessel functions ji(k q r) and spheri- 
cal harmonics Yz m (0, ip) and are characterized by the az- 
imuthal, polar and radial mode numbers m, t and q. For 
large m and £ the overlap factor in Eq. (|2j) breaks up into 
a product of the radial part cr r and the angular part ex- 
pressed in terms of the Clebsch-Gordan coefficients fl6j ] : 



cr = J ^ s ^i^l dV = cr r • {t 8 ,U',rn 8 ,mi\l p ,m p ). 



(3) 



The radial overlap factor <j r is found numerically, ap- 
proximating the large-index Bessel functions in vicin- 
ity of their q-th zero by Airy functions. It turns out 
that for large values of £ s ,i,p and m s ^ p accessible in our 
experiment, a r is nearly constant. Its dependence on 
q Sj i is shown in Fig. [T] for q p = 1. The maximum of 
a r « 230 cm -3 / 2 is reached for q s = qi = 1, however for 
other q s ^ the overlap remains significant. 
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FIG. 1: (Color online.) A) The radial overlap factor for the pump WGM nearest to the surface (q p = 1) and various signal 
and idler WGMs. B) The intensity distribution inside a WGM resonator for various mode numbers £, m q. C) The norma- 
lized angular overlap factor for an equatorial pump WGM and various signal and idler WGMs. 
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The Clebsch-Gordan coefficients arise when orbital 
momenta are added, e.g. in atomic physics. Their pres- 
ence in (j3]) indicates that in the PDC process in spherical 
geometry the pump photon's orbital momentum is shared 
by the signal and idler photons. The selection rules for 
this process replace the free-space phase matching con- 
ditions given by the second equation in (pQ): 

m s + uii = rn p 

\£ a -£i\<£ p <£a+£i (4) 

£ p + £ a + £i = 0, 2,4... 

To simplify the analysis of the orbital selection rules, 
we assume that in addition to being closest to the surface, 
which means that q p = 1 , the pump WGM is also equato- 
rial, i.e. £ p = m p . We evaluate the Clebsch-Gordan coef- 
ficients using Gaunt's formula [17[, treating large factori- 
als in Stirling's approximation. The coefficients coupling 
the equatorial signal and idler WGMs (£ s = ra s , £{ = mi) 
to the equatorial pump WGM (£ p = m p = m s H-ra^) only 
weakly depend on m p and are practically independent of 
how it is shared between the signal and idler photons: 

{m a ,mi\m s ,mi\m p ,m p ) « ^ (m p /7r 3 ) 1/4 . (5) 

In other words, the normalized angular overlap factor 
(m s + v S: rrii + Vi\ ra s , mi\m pj m p ) • 2 (tt 3 /rrip) 1 ^ 4 depends 
only on v s = £ s — m s and V{ = £ i — rrii 1 as shown in Fig. HJ 
Again, while the all-equatorial (i.e., v s = V{ = 0) WGMs 
coupling is the strongest, coupling of other signal and 
idler WGMs to the equatorial pump is also significant. 
Therefore we expect, and observe in experiment, several 
PDC channels for the same pump mode. Different cr and 
Q-factors of these modes lead to a different OPO thresh- 
old for each channel. For the most efficient all-equatorial 
case our analysis based on the measured Q s ,i ~ 3.4 x 10 7 
and Q p w 2.0 x 10 7 yields threshold power P th w 13 nW. 



For investigating these intriguing properties of a WGM 
OPO we built a setup schematically shown in Fig.[2j This 
setup is similar to the one used in [14[. It consists of a 
diamond turned and polished disk resonator made from 
a congruent 5% MgO-doped z-cut Lithium Niobate wafer 
of 0.5 mm height. The disk radius is R = 1.9 mm and 
the rim radius is 0.25 mm. 

A frequency-doubled Nd:YAG laser beam is coupled 
evanescently to our resonator via a diamond prism using 
frustrated total internal reflection. Adjusting the dis- 
tance between the prism and the resonator we change 
the overlap between their evanescent fields. This allows 
us to tune the coupling rate 70 of the WGM resonator, 
which would be equivalent to realizing mirrors with con- 
tinuously variable reflectivity in usual optical resonators. 

The optical axis of the crystal is chosen to coincide 
with the symmetry axis of the resonator. In this ge- 
ometry noncritical Type I PDC phase matching can be 
achieved between the extraordinary polarized pump (p) 
with a wavelength of 532 nm and ordinary polarized sig- 




FIG. 2: (Color online.) Experimental setup for PDC. 
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FIG. 3: (Color online.) WGM resonances at 532 nm and 1064 
nm wavelengths. 
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FIG. 4: (Color online, 
frequency detuning. 
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nal (s) and idler (i) with wavelengths near 1064 nm at 
a crystal temperature near 94° C. The triple resonance 
of the WGM resonator with the pump field and the two 
parametric fields is achieved by using the temperature 
disp ersion and the electro-optical effect, as discussed in 
|14j . Our measurement does not require any locking 
mechanisms, as the monolithic WGM resonators are re- 
markably stable. The non-degenerate signal and idler 
beams are coupled out by the same coupling prism, and 
separated for detection by a dispersion prism. 

To study the quality factor of the resonator, we mea- 
sured the critically coupled WGMs linewidths 7 at the 
pump and parametric wavelengths. In Fig. [3] we show 
two such modes with Q s ^ w 3.4 x 10 7 at 1064 nm, and 
Q p w 2.0 x 10 7 at 532 nm. The measurement of Q- 
factors reveals the resonator's intrinsic loss nature. The 
two dominant loss mechanisms in our resonator are mate- 
rial absorption and Rayleigh scattering at the resonator's 
defects and surface roughness. Radiation losses for this 
size of the resonator are negligible. Since the Rayleigh 
loss scales as the fourth power of the wavelength, one 
would expect a 16-fold increase of Q for twice the wave- 
length. Instead, we see an approximately 1.7- fold aug- 
mentation, attributed to stronger absorption of LiNbOa 
for shorter wavelengths [18|- Thus, our resonator quality 
factor is limited by material absorption only. 

The critital coupling contrast of 80% is limited by 
matching of the pump footprint spatial profile to the 
WGMs evanescent field in the coupling region. In the 
following we will operate with the "in-coupled" pump 
power, which is equal to 80% of the incident pump power. 

In contrast with Fig. [3j the pump WGM in Fig. [4] is 
phase-matched to the parametric WGMs. As soon as 
the in-coupled power exceeds the threshold value, the 
PDC fields are generated, which in turn affects the pump 
intensity. A series of such measurements for the same 
mode is summarized in Fig.[5j which shows the total OPO 
output power P s +i at the exact resonance as a function of 
the in-coupled pump power P. Note that the conversion 
efficiency is larger in the stronger coupled case. 

Following the theoretical analysis of the triply-resonant 



OPO dynamics jl9j |. we derive the following relation: 



Ps+i — &Pth 



TpoTo 
1p 7 



P/Pth 



(6) 



where j p o / "y p is the ratio of the pump WGM coupling rate 
to its linewidth, and 70/7 is a similar ratio for the signal 
and idler. Fitting the data in Fig. [5] with Eq. ([6]) we find 
Pth for an under coupled and a critically coupled pump. 
The lowest threshold power was estimated to be 6.7 fiW. 
While this is a significant improvement over previously 
reported values of 300 /xW for PDC Hj and 170 fiW for 
a four- wave- mixing process [2l|, the measured threshold 
is still much higher than expected. Partly this may be 
due to the discrepancy between the spherical approxima- 
tion adopted in our analysis and non-spherical geometry 
of our resonator. However, the main reason appears to 
be the non-optimal mode overlap, which can occur in 
many ways (see Fig. [I]). Therefore finding the optimal 
parametrically coupled WGMs should again increase the 
OPO efficiency by more than two orders of magnitude. 

For many spectroscopy and quantum optics applica- 
tions, narrow-line (ideally, single frequency mode) OPOs 
are desirable. We recorded the emission spectra of our 
WGM OPO for different phase matching conditions using 
an optical spectrum analyzer. Two such spectra labeled 
A and B are shown in Fig. [6] As seen from the figure, 
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FIG. 5: (Color online.) The measurement and theoretical 
fit for the resonant OPO output vs. in-coupled pump for a 
critically (inset) and weakly coupled resonator. 
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FIG. 6: (Color online.) Optical spectrum of the parametric 
light: two different mode pairs (A and B) are excited under 
different conditions. Inset: a zoom on one spectral component 
proves single-mode operation (solid and dashed lines represent 
diffirent conditions). 

the signal and idler are separated by more than 100 nm. 
Moreover, each parametric beam consists of a single fre- 
quency mode. To verify it, we recorded the signal and 
idler spectra with a resolution of approximately 2.7 GHz, 
which is sufficient to distinguish two WGMs separated 
by the signal and idler free spectral range (FSR) of ap- 
proximately 12 GHz. While recording this spectrum, the 
pump frequency sweeps across approximately 500 MHz, 
and the optical spectrum analyzer captures and retains 
all signal frequencies generated during the sweep. The 
result is shown in the inset of Fig. [6] The solid trace 
recorded over several minutes with no active locking of 
the resonator shows only one mode. For comparison, we 
have forced the resonator to "mode hop" across one FSR 
while recording the dashed trace, and verified that the 
two-mode operation would be clearly visible on the spec- 
trum analyzer. Thus we conclude that our OPO operates 
in a single mode regime. 

The single-mode operation is a consequence of the FSR 
dispersion A(A) = c/2itRn(\) which arises from the re- 
fractive index dispersion n(A). If the FSR dispersion is 



smaller than the signal and idler cavity bandwidth, i.e. if 
|A$ — A s | < 7, multi-mode PDC is allowed by the phase 
matching and energy conservation conditions. Then an 
FSR-spaced comb of pair-wise matched signal and idler 
frequencies is expected, centered at the degenerate PDC 
frequency. However, in a highly nondegenerate PDC the 
FSR dispersion is strong, |A$ — A s | > 7, and the single- 
mode regime is enforced by the energy conservation. 

The difference of the signal (A s ~ 1010 nm) and idler 
(Xi w 1120 nm) FSRs in our resonator can be found 
from the Sellmeier equation (22[: A^ — A s w 0.2 GHz, 
or from a direct measurement using the optical spectrum 
analyzer: A;— A s = 0.16±0.09 GHz. The theoretical and 
experimental results agree, and both greatly exceed even 
the overcoupled WGM's linewidth. Hence, the single- 
mode operation of our OPO is to be expected. 

To summarize, we have demonstrated a strongly non- 
degenerate, single mode, narrow line, high stability OPO 
with an extremely low threshold, based on a second-order 
nonlinear WGM resonator. Our analysis reveals the role 
of the photon orbital momentum conservation in the 
PDC phase matching in spherical geometry, explained 
the presence of multiple PDC channels with different ef- 
ficiencies, and predicted that even lower OPO thresholds 
could be achieved. This work illustrates the great poten- 
tial of WGM-based OPOs. They will allow for investigat- 
ing the rich dynamic phenomena in far above-threshold 
regime, and will provide the underlying technology for 
compact, stable, low-power tunable spectroscopy sources 
with narrow linewidth. Utilizing quasi phase matching in 
periodically poled resonators [TJ one could achieve 
flexibility in the selection of wavelengths. Finally, these 
resonators open the possibility to combine the optical 
nonlinear ity with optomechanics, where the crystalline 
materials are a promising candidate [24[. 
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